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The boson sampling problem has brought a lot of attention in the quantum information 
field because it is not efficiently solvable with a classical computer; nonetheless it can 
be implemented with linear optical interferometers with single-boson sources. Recently, 
we introduced a more general problem, the multi-boson correlation sampling problem , 
which allows us to take advantage of the multi-mode spectral distribution of the bosonic 
sources together with time-correlated measurements in order to achieve sampling not 
only over the output ports of the interferometer but also over the joint detection times. 
This problem was analyzed for both single-photon sources and thermal sources. In this 
work, we demonstrate that it is possible to use single qubit bosonic sources in order to 
sample not only over the described “space-time” degree of freedom but also over all the 
possible exponential number of multi-qubit output states. 

Keywords : quantum computing; boson sampling; multi-boson correlation sampling; 
multi-boson interference 


1. Introduction 

The boson sampling problem (BSP) 1,2 > 3 > 4 has recently raised a lot of attention in 
the research community leading to several proof-of-principle demonstrations 5 ’ 6,7,8 , 
as well as studies of its characterization 9 ’ 10,11 > 12,13,14,15,16 anf ] G f j^ s implemen¬ 
tation with different types of sources 17 : 18 . 19 > 2 °. 21 > 22 . 23 _ The BSP problem consists 
of sampling the probability distribution associated with N bosons at the output 
of a linear M-port interferometer, where the condition M N ensures that the 
probability for the bunching of bosons is reasonably small 1 . The BSP is hard to 
solve with a classical computer since it involves an exponential number of possible 
samples, each connected with the permanent of a random complex matrix. Indeed, 
the calculation of these permanents is a #P-problem 24 and therefore not efficiently 
computable with a classical computer. 

In recent works 25 i 26 . 2 c 28 ) we discussed the use of independent sources with 
arbitrary spectral distributions. Moreover, we analyzed the joint probability rate 
associated with the correlated detections of N bosons in a given sample of N output 
ports. Here, differently from the original formulation of the BSP, the sampling 
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problem relies on the times at which the bosons are detected. This leads to the more 
general multi-boson correlation sampling problem (MBCSP), which we analyzed in 
the case of Fock state sources 25 and thermal sources 26 . 

In this paper, we consider the case of arbitrarily polarized single-photon sources 
and we analyze the sampling problem over all the possible correlated single-photon 
detections both in time and polarization. By performing measurements in a given 
polarization base, it will be possible to sample over all the possible output JV-qubit 
states (corresponding to iV-binary-digit numbers) within each sample of N output 
ports. 

The paper is organized as follows: we give a brief summary of the MBCSP for 
single-photon sources in section 2; we consider the MBCSP based on polarization 
correlation measurements in section 3: after analyzing arbitrary polarization corre¬ 
lation measurements in section 3.1, we address the problem of the qubit sampling 
in section 3.2; we conclude with final remarks in section 4. 

2. The multi-boson correlation sampling problem (MBCSP) with 
multi-mode single-photon sources 

We summarize the MBCSP of a given order N with multi-mode single-photon 
sources as follows. First, we prepare a random linear interferometer (see Fig. 1) 
with M N ports. The interferometer input state |iF jv) is the 7V-photon state 
consisting of multi-mode single-photon states 29 

/ +oo 

dw^(w)alH|0), (1) 

-oo 

in each of N out of M given input ports s = 1 where £(ui) can be for 

simplicity a generic Gaussian narrow frequency distribution, and the vacuum in 
the remaining ports. Secondly, using time-resolving detectors, we sample from all 
possible correlated detection events in which N single photons are detected at only 
N of the M output ports at any possible sequence of correlated detection 

times for given integration intervals { At ^}. The total sampling time domain must 
be sufficiently large with respect to the coherence times of each single photon as 
well as the relative time delays of the photon pulses at the detectors. 

The rate of TV-fold joint detection events for a given sample d = 1,..., N of out¬ 
put ports is proportional to the TVth-order Glauber-Sudarshan correlation function 

30,31,32 

G W ({M) = (E^ ) (t 1 )...E«(t JV )Ei+ ) (t iV )...E^(t 1 )) (2) 

\ / I^jv) 

where E^ (t,i) denotes the positive/negative frequency parts of the field opera¬ 
tor E d(td) = E ^(td) + E ^\td) at the cith detector. These field operators are 
connected with the field operators at the input ports by a unitary operator U that 
describes the evolution in the linear interferometer. Equivalently, the interferometer 
can be described by a random unitary M X M matrix U^ MxM \ which we assume 
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Fig. 1. Implementation of the MBCSP of order N with a random linear interferometer with 
M N ports and N single-photon sources. After the evolution in the interferometer, described 
by a unitary random matrix iy{( MxM ) ? correlated detection events are recorded in N of the M 
output ports, labeled as d = 1,..., N. 


for simplicity to be frequency independent (any unitary matrix C an be 

implemented by using only passive linear optical elements 33 ). For a specific set of 
N occupied input ports and N output ports where a joint detection occurs, the 
N x N submatrix 


U (NxN) = 


u d . 


d= 1 , 
s= 1,. 


,,N' 

,N 


( 3 ) 


allows us to express the electric field operators at the detectors as linear combina¬ 
tions 


M 

E { d + \t d )=Y,Ud, s Ei+\t d ) (4) 

S= 1 

of the field operators e[ +) {td) at the sources. 

In Ref. 25 we found that the TVth-order Glauber-Sudarshan correlation function 
can be expressed in terms of matrices 

T (JVxJV) = [U d>a Xs(td))d=i,...,N, (5) 

s=l 

given by the Hadamard product of the interferometer matrix 7/( ArxAr ) in Eq. (3) 
and the time-dependent matrix [xs(td)]d=i,...,N, defined by the Fourier transforms 

a=l,N 

1 r + °° 

X s (t)=-7= 6(w) e lut duj (6) 

v "7T J — oo 


of the single-photon distributions £ s (u;), with s = 1 at times t = t d , with 
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d = 1, ...,7V . Indeed, in the narrow-bandwidth approximation Eq. (2) becomes 25 


G^\{t d })oc 


n [g<t(«),s x«(*<t(s))] 

5 = 1 


perT (7Vxiv) 


(7) 

( 8 ) 


Here, the permanent per 'f( NxN ') is defined as the sum over all the permutations a in 
the symmetric group Sn corresponding to all N\ amplitudes nfLi[Gcr(s),s X«(tr(s))] 
associated with each possible multi-photon quantum path from the N sources to 
the N detectors. 

In the case of ideal detectors we easily obtain, by using Eq. (8), the probability 


P (N) ({t d ,At d }) = 



f>t N +At N /2 

perT (JVxJV) 

/ dti .. 

. / dtN 

ltx-Ati/2 

J tj\f — At N /2 



(9) 


for an iV-folcl joint detection event in the output ports d = 1,2,..., N and in a given 
sample of time intervals {[t d — At d /2,t d + At d /2)} d=1 2 N . 


3. Boson Sampling based on polarization correlation 
measurements 

In this section, we show how we can take advantage of the polarization degree of 
freedom in order to implement a sampling problem not only in the output spatial 
channels and the detection times but also in the measured correlated polarizations 
of the output photons. 



Fig. 2. Implementation of the bosonic qubit sampling problem with a random linear interfe¬ 
rometer with M N ports and N single-qubit sources with polarization states |0 S ), where 
s = 1,2, ...,7V. After the evolution in the interferometer, described by a unitary random matrix 
Ij(MxM), correlated detection events both in time and polarization are recorded in N of the M 
output ports, labeled as d = 1,..., N. The polarization correlation measurements are performed 
by using polarization analyzers with generic orientations in front of the single photon detectors. 
Moreover, for a defined polarization base {| H), |V}}, it is possible for each port to measure a qubit 
either in the state | H) or |V) by using a polarization beam splitter and two detectors recording 
either the transmitted H-polarized qubit or the reflected V-polarized qubit. 
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For this reason we consider (see Fig. 2) N multi-mode single photon sources in 
the input ports s = 1, 2, ..TV, with polarizations 9 S1 corresponding, in a given base 
{\H), |V)} = ||0), 11)}, to N generic polarization qubits 

\0 8 ) = cos6» s |0) + sin(9 s |l), (10) 


leading to the iV-qubits state 


N 

Win) = 0 | 0 fl ) 

S=1 


21 ,... 


E 


r N 

q cos (6s 

- s =1 



( 11 ) 


with ||*«)• This state is given by the superposition of an expo¬ 
nential number 2 N of iV-qubit states |z) = with i = 0,1,2^ — 1 with 

amplitudes determined by the generic polarization angles 6 S , with s = 1,2, ...IV. 

After the evolution of these qubits in the interferometer, that we assume for 
simplicity to be polarization independent, polarization correlation measurements 
are performed by sampling over all M output ports. 


3.1. Correlation measurements at arbitrary polarization angles 


We consider first the case of correlation measurements at given angles <f> d of the 
polarization analyzers at each of the output ports, d = 1,2, ...,7V, where a joint- 
detection occurs. Here, it is useful to introduce the “effective” matrix 


gOd A) 


cos {(j) d - Q s )U d , s 


d=l,...,N 

s=l,...,N 


( 12 ) 


obtained by the Hadamard product of the matrix lf ( NxN ) i n Eq. (3) and the polar¬ 
ization matrix [cos(</f>d — 0 s )]d=i,...,iV) defined by the source polarizations {0 S } and 

s=i ,’jv 

the detected polarizations {4> d }. 

By substituting the polarization-independent interferometric matrix in 

Eq. (3) with the “effective” matrix (12) the TV-order correlation function in Eq. (7) 
becomes 


G (Ar) ({f d ,<M;{0 s }) oc 


per \13 


{4>dfis) 

d,s 


Xs(td)] 


d=l,... t N 

s=l,...,N 


2 


In an analogous way, Eq. (9) corresponds to 


p (N) ({t d ,At d ,M-,m) = 


nti+Ati /2 


P er [l3 i J , s’ es) Xs(td)\d= 1.JV 

/ dt i .. 

. / dtpf 

J ti — Ati/2 

1 

to 

s=l 


(13) 


(14) 


representing the TV-fold joint detection probability of measuring the N qubits with 
polarization angles <f> d , for a given sample d = 1,2, ...N of output modes and detec¬ 
tion time intervals {[t d — At d /2,t d + At d /2]}. 

In conclusion, we have obtained a compact expression in terms of permanents for 
the probability of detecting the N qubits in each possible sample of output ports, 
















FILE 


June 17, 2015 0:26 WSPC/INSTRUCTION 

QubitBosonSampling_tamma 


6 Vincenzo Tamma 

depending not only on the detector integration intervals but also on the general 
orientation of the polarization analyzers in front of the detectors. 


3.2. Sampling of photonic qubits 

In this section, we finally describe how the quantum network in Fig. 2 can be used 
to implement a multi-boson correlation sampling problem, where the sampling oc¬ 
curs not only over the output spatial channels and the joint-detection times but 
also over all possible output TV-photon polarization outcomes, with respect to the 
base {| H), |V)} = {|0), 11)}. In order to perform correlated polarization measure¬ 
ments in such a base we use, for each output port, a polarization beam splitter 
and two detectors recording either a transmitted TJ-polarized qubit or a reflected 
H-polarized qubit. Thereby, it is possible for each port to measure a qubit either 
in the state |0) or |1). Indeed, our device is sampling over all the possible qubit 
correlation measurements with outcomes {|jd)}d=i,...,jv, where j d = 0,1. 

Depending on the measured qubit states {|jd)}d=i,...,jv for an y possible “spatial” 
sample d = 1,2,..., N the “effective” matrices (12) reduce to 


fiiidfis) 


cos (-jd - 0 s )Ud,s 


d= 1,... 
s= 1,... 


.,JV» 
,N 


(15) 


and the polarization correlation probability in Eq. (14) becomes 


P (N) ({t d ,At d ,j d };{6 s }) = 



nt n ~\~ At n / 2 

per [B^ d / s) x s (t d )\ d = 

/ dU •• 

. / dt N 

/ti-Ziti/2 

o tjy — At tv /2 

’ s=l,...,N 


for the IV-qubit input state in Eq. (11). Here, the TV! amplitudes in the permanent 
in Eq. (16), associated with all the interfering multi-photon quantum paths, depend 
not only on the interferometer evolution, the source frequency distribution and de¬ 
tection time intervals, but also on which of the possible 2 N qubit states \j\---jd---3N) 
is detected in the qubit sampling problem. 


4. Conclusion 

In this paper, we generalized the multi-boson correlation sampling problem to pho¬ 
tonic qubits, by taking advantage of the polarization degree of freedom of the input 
photons. Given a linear interferometer with single-photon qubit sources, we defined 
a more general problem based on the sampling of the output probability distribu¬ 
tion not only in the spatial channels and associated detection intervals but also in 
all possible qubit measurement outcomes, exponential in number, in a given polar¬ 
ization base. We found, for each possible degree of freedom in the sampling, i.e. 
spatial channels, time, and polarization, a corresponding matrix. It is the modulo 
squared permanent of the Hadamard product of these matrices which determines 












FILE 


June 17, 2015 0:26 WSPC/INSTRUCTION 

QubitBosonSampling_tamma 


Sampling of bosonic qubits 7 


the detection probability rate for each possible sample. Finally, the probability asso¬ 
ciated with each correlation measurement emerges naturally from integrating over 
the detection time intervals associated with the time sampling. 

The result obtained so far can be easily extended to thermal sources 26 and, 
in general, from photonic networks to atomic interferometers with bosonic qubits, 
where, for example, the spin degree of freedom takes the place of the polarization. 
Although for simplicity we have addressed the case of polarization-independent 
evolution through the interferometer, our work can be generalized to any non-trivial 
interferometric evolution in the multi-qubit Hilbert space. 

This research may pave the way towards the developments of quantum networks 
able to achieve quantum speed-up in useful computational problems beyond boson 
sampling. 
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